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In  this  study,  we  investigate  the  air- water  two-phase  flow  in  a  single  flow  channel  of  polymer  electrolyte 
membrane  (PEM)  fuel  cells.  In  the  ex  situ  study,  both  straight  and  serpentine  channels  with  various  gas 
diffusion  layer  (GDL)  surfaces  are  studied.  Focus  is  placed  on  the  two-phase  flow  patterns,  which  are 
optically  characterized  using  a  microscope  with  a  high-resolution  camera,  and  the  two-phase  pressure 
amplifiers.  We  find  that  the  GDL  surface  properties  slightly  affect  the  flow  pattern  and  two-phase  pressure 
amplifier  in  the  flow  field  configuration.  Flow  pattern  transition  occurs  at  the  superficial  gas  velocity  of 
around  1ms-1,  and  the  pressure  amplifier  can  reach  as  high  as  10.  A  two-fluid  model  is  also  presented 
together  with  one  dimensional  (1-D)  analytical  solution,  and  acceptable  agreement  is  achieved  between 
the  model  prediction  and  experimental  data  at  high  gas  flow  rates. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cells,  which  convert 
the  chemical  energy  stored  in  hydrogen  fuel  directly  and  efficiently 
to  electrical  energy  with  water  as  the  only  by-product,  have  the 
potential  to  reduce  our  energy  use,  pollutant  emissions,  and  depen¬ 
dence  on  fossil  fuels  [1,2].  Gas  flow  channels  are  an  important  fuel 
cell  component  for  reactant  supply  and  water  removal.  The  surface 
of  gas  diffusion  layers  (GDL)  constitutes  one  wall  of  a  gas  flow  chan¬ 
nel,  playing  an  important  role  in  the  water  management  of  fuel  cells 
[3,4]. 

Air-water  two-phase  flow  has  been  the  focus  of  many  stud¬ 
ies.  An  important  issue  is  to  classify  different  types  of  flows.  Flow 
regimes  can  be  displayed  in  a  flow  regime  map  [5-7]  with  superfi¬ 
cial  gas  and  liquid  velocities  as  coordinates.  Typical  patterns  include 
bubbly,  plug,  slug,  slug-annular,  dispersed  and  annular  flows.  The 
frictional  pressure  drop  APf  is  another  important  parameter,  which 
can  be  characterized  using  the  two-phase  flow  multiplier  defined 
as  the  ratio  of  the  two-phase  pressure  drop  to  the  single-phase  one 
[8-10].  In  PEM  fuel  cells,  the  presence  of  liquid  water  increases  A Pf 
in  gas  channels  [11-15].  Wang  et  al.  explained  that  the  increased 
A  Pf  due  to  water  accumulation  is  one  important  cause  for  flow 
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mal-distribution,  and  that  the  two-phase  flow  multiplier  may  be  a 
reliable  indicator  of  water  buildup  in  channels  [14,15]. 

Liquid  water  in  channels  can  be  visualized  using  optical  methods 
[16-22].  Tuber  et  al.  were  among  the  first  groups  to  analyze  water 
accumulation  in  a  gas  channel  using  a  transparent  fuel  cell  [16].  Lu 
et  al.  designed  an  experimental  setup  with  0.7  mm  inner  diameter 
channels  and  showed  three  flow  regimes  [19].  Neutron  radiogra¬ 
phy  investigations  reported  that  higher  water  retention  may  be 
found  at  the  U-bends  of  the  serpentine  channels  [23],  and  that 
the  use  of  hydrophobic  coating  on  the  gas  channels  leads  to  larger 
water  retention.  Modeling  has  been  conducted  by  many  studies 
in  the  area  of  two-phase  channel  flows.  Most  used  the  volume  of 
fluid  (VOF)  method,  which  tracks  the  liquid-air  interface.  However, 
the  VOF  approaches  are  usually  numerically  expensive,  therefore 
it  is  difficult  to  extend  it  to  couple  with  the  electrochemical  kinet¬ 
ics  and  other  transport  mechanisms  in  a  fuel  cell  model.  In  this 
regard,  an  engineering  sub-model  that  is  able  to  describe  the  impor¬ 
tant  features  of  the  two-phase  channel  flows  and  can  be  efficiently 
incorporated  into  a  fuel  cell  model  is  highly  desirable. 

In  this  work,  ex  situ  experiment  is  conducted  to  examine  the 
characteristics  of  the  two-phase  flows  in  a  single  flow  channel  of  a 
fuel  cell.  A  two-fluid  model  is  presented  and  its  analytical  solution 
to  the  experimental  case  is  derived.  The  model  prediction  com¬ 
pares  with  the  experimental  data  in  terms  of  two-phase  pressure 
amplifier.  In  addition,  though  many  in  situ  studies  were  conducted 
previously  to  reveal  the  flows  inside  fuel  cell  channels,  it  is  dif¬ 
ficult  to  characterize  the  two-phase  flows  for  modeling  purpose 
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Fig.  1.  (a)  Schematic  of  the  experimental  setup;  (b)  detailed  view  of  the  experimental  device;  (c)  the  straight  channel  flow  field  (only  the  3rd  channel  will  be  tested);  and  (d) 
the  serpentine  flow  field. 


because  the  in  situ  flow  rate  varies  along  the  channel  due  to  local 
water  production  by  fuel  cells.  Ex  situ  experimental  data  under  well 
controlled  flow  conditions  are  particularly  important  for  the  initial 
stage  of  model  development. 

2.  Two-phase  flows  in  the  flow  fields  of  PEM  fuel  cells 

In  fuel  cells,  water  is  produced  in  the  cathode.  Liquid  emerges 
when  the  partial  vapor  pressure  reaches  the  saturated  one,  and 
eventually  enters  the  channel  via  the  GDL,  resulting  in  air-water 
two-phase  channel  flow.  The  water  in  channels  can  accumulate  to 
a  degree  that  greatly  affects  reactant  supply.  The  water  generated 
rate  can  be  expressed  using  the  Faraday  law: 

Sw=T  (!) 


be  developed  using  the  concept  of  the  two-phase  flow  multiplier 
4>l  [8,9]: 


-(dP/dz)fr2<p 

0C  -O  P/te)fic 


(5) 


Chisholm  examined  the  Lockhart-Martinelli  correlating  proce¬ 
dure  and  gave  the  pressure  multiplier  for  pipe  two-phase  flows 
[10]: 


4>l  =  l  +cx+x2 


(6) 


where  X  is  the  Lockhart-Martinelli  parameter.  Mishima  et  al.  indi¬ 
cated  that  the  coefficient  C  can  depend  on  the  micro-channel 
diameter  [24]: 

C  =  21(l  +e"0319Dh)  (7) 


The  net  water  transport  coefficient  a  is  frequently  adopted  to 
calculate  the  net  water  gain  in  the  cathode: 

Scw  =  (  l+2a)T  (2) 

Assuming  all  the  water  addition  is  in  liquid  form  when  enter¬ 
ing  channels,  integrating  the  above  flux  over  a  fuel  cell  yields  the 
liquid  flow  rate  in  the  cathode  channel.  The  inlet  gas  flow  rate  is 
determined  by  the  stoichiometric  ratio  §c- 


I  inlet 


£cfAn 

4FCq2Ac 


(3) 


where  Ac  and  Am  are  the  flow  cross-sectional  areas  of  the  cathode 
gas  channels  and  the  membrane  area,  respectively.  For  a  typical 
condition  of  2  atm,  80  °C,  1.5  stoichiometry,  and  1 A  cm-2,  the  gas 
velocity  is  in  the  magnitude  of  1ms-1. 

In  gas  channels,  flow  regimes  can  be  displayed  using  the  super¬ 
ficial  gas  ( UG )  and  liquid  ( Ul )  velocities  as  coordinates: 


UG  =  and  UL 

Pg 


mxL 

Pl 


(4) 


where  xG  and  xL  are  the  qualities  of  gas  and  liquid,  respectively,  and 
fn  the  mass  flux.  The  frictional  pressure  drop  (A Pf)  is  an  important 
parameter  characterizing  two-phase  flows,  and  its  correlation  can 


where  the  hydraulic  diameter  is  defined  as: 

4  x  cross  sectional  area 

Dh  =  - : - : - : -  (8) 

channel  perimeter 

3.  Experimental 

In  experiment,  the  air-water  two-phase  flow  was  generated 
in  a  horizontal  channel  with  controlled  flow  rates,  see  Fig.  1.  The 
cross  section  of  flow  channels  is  rectangular  with  a  1.6  mm  width 
and  1  mm  height.  The  distance  between  the  air  inlet  and  outlet 
is  50  mm  for  the  straight  channel.  The  header  has  a  dimension  of 
~10mm  x  5  mm  x  1  mm.  Three  types  of  GDL  surfaces  are  consid¬ 
ered:  smooth  aluminum,  carbon  cloth,  and  carbon  paper.  Both  the 
carbon  cloth  and  paper  are  treated  hydrophobic  with  30%  PTFE 
content.  The  side  walls  of  the  channels  are  made  of  aluminum, 
while  the  top  uses  a  hydrophilic  polycarbonate  plate  to  facilitate 
optical  access.  Fig.  2  displays  the  water  droplet  shapes  over  the 
aluminum  (contact  angle  ~50°)  and  polycarbonate  plates  (contact 
angle  ~70° ),  respectively.  Fig.  1  (b)  shows  a  detailed  view  of  the  test 
section.  The  gaseous  phase  was  fed  by  a  mass  flow  controller.  Water 
was  introduced  to  the  channel  through  a  ~0.1  mm  diameter  circular 
hole  located  on  the  bottom  of  the  channel  (GDL  side),  close  to  the  air 
inlet.  Water  flow  rate  was  controlled  by  a  Chemyx  Inc.  Fusion  400 
micro  pump.  Absolute  and  relative  pressure  measurements  were 
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Fig.  2.  The  pictures  of  droplets  over  polycarbonate  (left)  and  aluminum  (right)  surfaces. 


taken  by  a  pressure  transducer.  An  overhead  microscope  and  cam¬ 
era  were  used  to  record  the  images  of  the  flow  regimes.  Fig.  1(a) 
shows  the  schematic  of  the  overall  experimental  setup.  The  super¬ 
ficial  gas  velocity  in  the  flow  channel  Uq  varies  from  0.44  m  s-1  to 
5.56  ms-1,  and  the  liquid  one  Ui  from  l.lxl  0-4  to  5.5  x  1 0-3  m  s-1 
for  most  cases.  It  is  worthy  to  mention  that  the  considered  region 
includes  a  flow  channel  and  headers  (also  called  manifolds)  for  the 
pressure  drop  measurement.  Further  in  real  PEFC  operation,  liquid 
water  comes  out  of  the  GDL  surface  via  multiple  pores.  We  hope 
through  this  single-pore  case,  i.e.  the  single  liquid  injection  port, 
we  can  verify  the  model  and  further  extend  it  to  our  future  study 
on  more  realistic  condition. 


Within  the  channel,  water  transport  takes  place  in  either  gas  or 
liquid  phase.  Note  that  the  liquid  phase  is  purely  water.  One  can 
derive  the  water  conservation  as: 

V  .  (ugCgH2°  +  u,^)  =  V  .  (£SgDg2°VCg2°)  (15) 

For  analysis  purpose,  we  consider  one  dimension  ( 1  -D)  along  the 
channel  or  x  direction.  We  define  a  dimensionless  x*,  scaled  by  the 
length  of  the  channel,  as  the  location  the  liquid  is  injected.  Before 
x*,  the  air  flow  is  in  the  single-phase  flow  regime.  After  x*,  liquid 
starts  to  form,  resulting  in  two-phase  transport.  One  can  integrate 
both  mass  conservation  equations  from  x*  to  x 


4.  A  two-fluid  model 

Due  to  the  mini/micro-scale  nature  of  the  channel  dimension, 
the  internal  two-phase  flow  can  be  treated  as  flows  in  porous  media 
[14,15].  The  mass  exchange  occurs  at  the  phase  interface  through 
water  evaporation  or  condensation.  The  mass  conservation  of  each 
phase  can  be  written  as  follows  [15]: 

Gas  phase  mass  conservation  :  V  •  ( pgug )  =  (9) 

Liquid  phase  mass  conservation  :  V  •  (p/U/)  =  -S^  (10) 


i  px 

PgUg  =  Pg(X*)Ug(x*)  +  ^  J  (Sm  +  psm)dx  (16) 

PiUi  =  L  f  (-s£  +  (l  - P)Sm)dx  (17) 

Lz  Jx* 

where  Sm  and  are  the  mass  addition  and  the  mass  fraction  of  water 
addition  in  form  of  vapor  from  the  channel  wall,  respectively.  Note 
that  in  this  case  is  zero  in  the  entire  region  except  the  location 
of  x*.  By  neglecting  diffusion,  applying  the  same  integration  over 
the  water  conservation  equation  yields: 


where  ug  and  ut  are  the  superficial  velocities  of  gas  and  liquid, 
respectively.  The  gas  density  can  be  calculated  through  the  con¬ 
stituent  species: 

pg  =  ^JWkQ  (11) 

k 

Here,  k  represents  species.  In  the  cathode,  there  are  nitrogen,  water, 
and  oxygen. 

Flow  in  the  porous  media  is  driven  by  the  pressure  gradient  and 
body  forces,  balanced  with  the  shear  stress  at  the  solid  surface.  By 
adopting  the  Darcy’s  law  and  neglecting  the  gravitational  force,  the 
momentum  equations  can  be  written  as: 

k  K 

Gas  phase  momentum  conservation  :  pgug  =  ~^—VPg  (12) 

vg 

k  iK 

Liquid  phase  momentum  conservation  :  p/U/  =  VP/  (13) 

vi 

The  relative  permeabilities,  krg  and  krh  are  defined  as  the  ratio 
of  the  intrinsic  permeability  of  liquid/gas  phase  to  the  total  per¬ 
meability  of  a  porous  medium.  Physically,  it  describes  the  extent 
to  which  one  fluid  is  hindered  by  others  in  pore  spaces  and  can  be 
approximated  as  a  function  of  liquid  saturation.  A  set  of  functions 
have  been  frequently  used  in  the  fuel  cell  literature  for  the  relative 
permeabilities: 

kri  =  s"k  and  krg  =  (1  -se)nk 

where  the  effective  saturation  se  =  ^ — —  (14) 

1  —  sir 


UsC”2° + u« =  %(**)cgH§a°t  +  rz[  sli2°d*  (18) 

where  SH2°  is  the  water  addition  rate  from  the  channel  wall.  Note 
that  in  this  case  SH2°  is  zero  in  the  entire  region  except  the  location 
ofx*. 

Assuming  the  interfacial  area  between  liquid  and  gas  phases  is 
large,  the  equilibrium  between  liquid  and  gas  phases  holds  true,  i.e. 
C|*2°  equals  to  C^2°t.  Assuming  a  constant  gas  density  pgl  one  can 
solve  the  phase  superficial  velocities 

=  LzPl{C"fa(M^o  .  pg)  (Wt  J_t  Smdx  ~  Ps  £  SH2°d*)  (19) 
and 

Ur  =  Ugi**]  "  Lz(CH20tMH2o_Pg)  ( l  -  MH2°  £  SH20£K) 

(20) 


Further,  the  Darcy’s  laws  for  both  phases  can  be  written  in  one 
dimension  as 


Ua  =  ~ 


krgK  dPg 
fig  dx 


u  =  _l<rl^  dPi  _  j  urg 

1  fii  dx  fii  l  dx 


krtK  (  dPg  dPc 
dx 


(21) 
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Fig.  3.  Flow  patterns  of  two-phase  flow  in  the  straight  channel  without  GDLs  (a).  From  the  above  to  the  bottom,  the  flow  patterns  are  annular,  wavy-annular,  wavy,  and 
slug-annular  flows  (b). 


Combining  these  two  equations  to  eliminate  (dPg/dx)  and  then 
substituting  the  Leverett  J-function  yields: 


^sKiu  |  /cri(/fe)1/2gcos(ec)  dj  ds, 
jjji  krg  fiiLx  ds  dx 


(23) 


The  above  equation  can  be  solved  numerically.  To  simplify  the 
analytical  solution,  one  can  assume  the  liquid  transport  along  the 
channel  is  primarily  driven  by  the  two-phase  interaction  through 
the  shear  stress,  and  the  capillary  action  is  small  and  negligible  [15]. 
Therefore  the  above  equation  can  be  resolved  explicitly: 


1  +  gUgP/'nk'1 


(24) 


Assuming  a  constant  air  superficial  velocity,  the  pressure  scale 
(i.e.  the  denominator  in  the  middle  term  of  the  above  equation) 
becomes  the  pressure  drop  for  a  single-phase  flow.  The  above  equa¬ 
tion  can  further  be  rearranged  as 


(28) 


It  is  worth  to  note  that  the  above  dimensionless  pressure  is 
defined  in  the  same  way  as  the  two-phase  pressure  amplifier 
in  Eq.  (5).  Detail  regarding  the  two-fluid  flow  model  and  analysis 
can  be  found  in  Ref.  [15]  for  fuel  cells. 


5.  Results  and  discussion 


Once  the  saturation  profile  is  available,  the  relative  permeabili¬ 
ties  can  be  calculated,  which  are  further  substituted  to  the  pressure 
equations 


/  llgUgdx 

L  w 

(25) 

r x 

(26) 

Pg  -  Pg(x*)  -  Lx 
P,=Ps(x*)-Lx 


The  total  pressure  drop  of  the  gas-phase  along  the  porous  chan¬ 
nel  will  be: 


APg  =  , 

(MgU. 


a* _ =  r^dk+  f  ^  a 

g.inC  )/K  J0  Ugjn  J-t  Ug  ln  kn 


- — dx 
rg 


(27) 


Fig.  3(a)  displays  the  two-phase  flow  regime  map  for  the  straight 
channel  without  GDLs,  i.e.  the  GDL  side  is  an  Aluminum  surface. 
Fig.  3(b)  shows  several  typical  flow  patterns  observed.  It  can  be  seen 
that  with  decreasing  gas  flow  rate  the  two-phase  flow  experiences 
annulus,  wavy  annulus,  wavy,  and  slug  in  the  range  of  considered 
operating  condition.  The  shear  stress  exerted  by  the  gas  stream  is  a 
major  driven  force  for  liquid  flow.  When  there  lacks  of  a  sufficient 
gas  flow,  liquid  accumulates,  forming  slug  flow.  When  the  gas  flow 
becomes  strong,  the  two-phase  stream  is  stable  annulus  with  liq¬ 
uid  over  the  channel  walls.  Fig.  4  presents  comparison  with  other 
observations  [19,25],  showing  an  acceptable  agreement. 

Fig.  5  displays  the  two-phase  pressure  amplifier  profiles  and 
comparison  with  the  Lockhart-Martinelli  correlation  and  model 
prediction.  The  pressure  amplifier  represents  the  effect  of  liquid 
presence  on  the  pumping  pressure  required  to  drive  the  gas  flow. 


Fig.  4.  Comparison  with  the  results  obtained  by:  (a)  Lu  et  al.  [19],  using  a  0.7  mm  inner  diameter  circular  tube;  and  (b)  Yang  et  al.  [25],  using  a  1  mm  inner  diameter  circular 
tube. 
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Fig.  5.  0^  as  a  function  of  Ug  and  UL  from  model  prediction  and  experiment  (a);  and  comparison  between  the  experimental  and  Lockhart-Martinelli  values  of  0^  (b). 


Fig.  6.  Flow  patterns  in  the  straight  channel  with  carbon  cloth  (a)  and  paper  (b). 


It  can  be  seen  from  Fig.  5(a)  that  the  value  of  the  amplifier  factor  is 
close  to  1  when  the  gas  flow  rate  is  high,  and  increases  to  around 
10  at  low  gas  flow  rates.  The  liquid  flow  rate  also  affects  the  ampli¬ 
fier.  The  model  predictions  are  also  provided  for  higher  gas  flow 
rates.  The  model  parameters  are  given  in  Table  1.  It  can  be  seen 
that  the  model  predictions  are  in  line  with  the  experimental  data 
for  the  considered  region.  For  lower  gas  flow  regions,  the  model 


parameters  must  be  changed,  for  different  liquid  flow  rates  and 
GDL  surfaces,  in  order  to  fit  with  the  experimental  data,  which  is 
likely  due  to  the  change  in  two-phase  flow  patterns.  Further  efforts 
are  required  to  examine  the  various  model  parameters.  We  will 
present  the  results  in  our  future  work. 

Fig.  5(b)  compares  the  two-phase  amplifier  with  the 
Lockhart-Martinelli  correlation,  indicating  a  reasonably  good 


Fig.  7.  Two-phase  pressure  amplifiers  for  straight  channels  with  carbon  cloth  (a)  and  paper  (b). 
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Fig.  8.  Two-phase  flow  dynamics  near  the  channel-header  connection  [30]. 


Table  1 

Geometrical,  physical  and  model  parameters. 


Quantity 

Value 

Gas  channel  depth/width 

1.0/1.6mm 

Viscosity  of  liquid  water,  u/ 

0.000894  kg  m-1  s”1 

Viscosity  of  air,  iig 

1.86  x  10-6  kgirr1  s_1 

Residual  liquid  saturation,  sir 

0 

Exponent  in  the  relative  permeability,  n/< 

2 

agreement.  We  used  the  Mishima  and  Hibild  approach  in  order 
to  calculate  the  coefficient  C  in  the  Chisholm  equation  [26].  It  is 
worthy  to  mention  that  the  gas  pressure  drop  is  closely  related  to 
the  pumping  power  required  to  feed  the  reactants  [27].  Another 
issue  closely  related  to  pressure  drop  is  flow  mal-distribution 
phenomena  occurring  among  parallel  channels.  Assuming  two 
channels  in  parallel  follow  the  same  trend  of  two-phase  pressure 
amplifier,  the  two  channels  will  have  totally  different  air  flow  rates 


when  subjecting  to  different  flow  regimes.  Other  investigations 
[28,29]  also  presented  similar  results  of  the  two-phase  pressure 
amplifier. 

Figs.  6  and  7  display  the  two-phase  flow  regimes  and  pressure 
amplifiers  for  carbon  cloth  and  paper.  When  using  a  hydrophobic 
surface  in  the  GDL  side,  we  found  that  visually  distinguishing  the 
wavy,  annulus,  and  wavy  annulus  patterns  becomes  difficult  due  to 
the  distinct  wall  properties  between  the  GDL  surface  and  channel 
side  walls:  the  liquid  free  surface  observed  through  the  top  trans¬ 
parent  plate  can  be  disturbed,  i.e.  wavy,  even  the  flow  still  resides  in 
the  stable  annulus  regime.  Therefore,  for  these  two  cases  we  only 
distinguish  the  slug-annulus  and  wavy-annulus  patterns,  i.e.  the 
wavy-annulus  pattern  lumps  the  ones  of  wavy,  wavy-annulus,  and 
possible  annulus.  It  can  be  seen  that  the  two-phase  flow  regimes 
are  similar  between  the  carbon  paper  and  cloth.  The  flow  transition 
takes  place  at  UG  of  around  1ms-1.  The  model  predictions  are  also 
plotted  for  comparison.  It  can  be  seen  the  agreement  is  acceptable. 


Fig.  9.  Two-phase  flow  visualizations  for  the  smooth  aluminum  surface:  (a)  and  (b);  for  the  carbon  paper  surface  (c);  and  for  the  carbon  cloth  surface  (d). 
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Fig.  10.  Two-phase  flow  patterns  and  pressure  amplifiers  for  a  serpentine  channel  with  the  carbon  cloth  (a)  and  (b),  and  carbon  paper  (c)  and  (d),  respectively. 


Note  that  the  model  does  not  distinguish  the  GDL  surface  proper¬ 
ties,  therefore  the  model  prediction  results  are  the  same  as  that  in 
Fig.  5. 

In  comparison  of  the  three  cases,  it  can  be  seen  that  the  differ¬ 
ence  is  small,  and  the  difference  occurs  primarily  in  the  transition 
region.  This  is  likely  due  to  the  difference  in  the  surface  properties, 
either  wettability  or  roughness.  In  addition,  we  also  find  that  at  the 
connection  between  the  flow  channel  and  header,  liquid  accumu¬ 
lates  in  the  header  side  and  can  block  the  channel  at  low  gas  flow 
rates,  see  Fig.  8,  which  shows  the  liquid  surface  dynamics  near  the 
channel  outlet.  Wang  et  al.  indicated  liquid  will  buildup  near  the 
interface  between  small  and  large  channels  due  to  heterogeneity 
[14].  This  liquid  build  up  may  block  the  small  channel  due  to  the 
capillary  action.  Note  that  many  designs  of  optical  fuel  cells  have 
similar  channel-header  configuration.  As  to  the  two-phase  pres¬ 
sure  amplifier,  the  profiles  exhibit  similar  trend  among  the  three 
cases. 

Another  popular  flow  field  in  PEM  fuel  cells  is  the  serpentine 
pattern,  see  Fig.  1(d).  Liquid  water  may  be  trapped  at  the  corners  of 
the  U-turns,  see  Fig.  9(a)  and  (c).  The  trapped  water  can  be  stripped 
off  by  gas  flow,  influencing  the  downstream  flow,  see  Fig.  9(b). 
For  comparison  purpose,  we  put  together  the  cases  for  the  carbon 
paper  and  cloth,  see  Fig.  1 0.  Again  similar  to  the  single-channel  case 
the  flow  transition  occurs  around  UG  =  1 .0  m  s-1 ,  indicating  that  the 
liquid  accumulate  near  the  corner  has  a  small  impact  on  the  two- 
phase  flow  patterns,  and  the  two-phase  pressure  amplifiers.  This 
is  possibly  due  to  the  fact  that  liquid  accumulation  at  the  U-turn  is 
relatively  slow  given  the  small  liquid  water  injection  rate  (usually 
the  case  in  a  fuel  cell)  and  the  stripped-off  liquid,  if  forming  slug, 
can  be  rapidly  removed  out  of  the  short  channel  by  gas  flow  in  the 
timescale  around  0.01-0.1  s,  therefore  imposing  a  small  influence 
on  the  averaged  two-phase  pressure  amplifiers.  The  model  predic¬ 


tions  are  also  plotted  in  Fig.  10(b)  and  (d).  Again  the  agreement  is 
acceptable. 


6.  Conclusions 

This  study  conducted  ex  situ  investigation  and  model  analysis  of 
air-water  two-phase  flows  in  a  single-channel  flow  field  of  PEM  fuel 
cells.  The  experiment  examined  two-phase  flow  patterns  and  pres¬ 
sures.  We  found  altering  the  GDL  surface  properties  may  slightly 
change  the  flow  patterns,  and  the  liquid  water  can  accumulate  near 
the  header-channel  interface,  affecting  the  channel  flow.  Flow  pat¬ 
tern  transition  occurs  around  the  gas  velocity  of  1  ms-1  for  the 
considered  channels,  and  the  two-phase  pressure  amplifiers  can 
reach  about  1 0  at  the  gas  velocity  of  0.5  ms-1,  and  are  reduced  to  1 
when  the  gas  flow  rate  is  sufficiently  high.  In  the  serpentine  config¬ 
uration,  we  found  that  the  water  trapped  at  the  corner  of  the  U-turn 
has  a  small  impact  on  two-phase  flow  pattern  and  pressure  ampli¬ 
fier  for  a  single  channel  flow.  A  two-fluid  model  was  presented  for 
the  experimental  case  with  analytical  solution  derived.  By  compar¬ 
ison,  we  obtained  an  acceptable  agreement  with  experimental  data 
at  high  gas  flow  rates. 
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